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The structure of strange dwarfs and that of hybrid stars with same baryonic number is 
compared. There is a critical mass (M ^ 0.24Mq) in the strange dwarf branch, below 
which configurations with the same baryonic number in the hybrid star branch are more 
stable. If a transition occurs between both branches, the collapse releases an energy of 
about of 3 X 10^'' erg, mostly under the form of neutrinos resulting from the conversion of 
hadronic matter onto strange quark matter. Only a fraction (~ 4%) is required to expel 
the outer neutron-rich layers. These events may contribute significantly to the chemical 
yield of nuclides with yl > 80 in the Galaxy, if their frequency is of about one per 1500 
years. 
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1. Introduction 

It has been suggested already in the 70's ^'^ that, because of the extreme densities 
reached in the core of neutron stars (NS), hadrons can melt, creating a deconfined 
state dubbed the "quark-gluon plasma". Stars with a deconfined core surrounded 
by hadronic matter are called hybrid stars (HS), whereas objects constituted by 
absolutely stable strange quark matter are christened "strange stars" (SS) 

A pure SS is expected to have a sharp edge with a typical scale defined by the 
range of the strong interaction. It was pointed by Ref. 4 that surface electrons may 
neutralize the positive charge of strange quark matter, generating a high voltage 
dipole layer with an extension of several hundred fermis. As a consequence, SS are 
able to support a crust of nuclear material, since the separation gap created by the 
strong electric field prevents the conversion of nuclear matter into quark matter. 
The maximum density of the nuclear crust is essentially limited by the onset of 
the neutron drip {pd « 4 x 10^^ gcm^^), since above such a value free neutrons 
fall into the core and are converted into quark matter. Stars with a strange matter 
core and an outer layer of nuclear material, with dimensions typically of white 
dwarfs, are dubbed "strange dwarfs" (SD) The stability of these objects was 
examined in Ref. 7, where the pulsation frequencies for the two lowest modes n ~ 
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0, 1, as a function of the central density, were studied. More recently, the formation 
of deconfined cores has been considered in different astrophysical scenarios. The 

hadron-quark phase transition induces a mini-collapse of the NS and the subsequent 
core bounce was already invoked as a possible model for gamma-ray bursts (GRB). A 
detailed analysis of the bounce energetics ^ has shown that the relativistic (F > 40) 
fraction of the ejectum is less than 10^^ org, insufficient to explain GRBs. However, 
as the authors have emphasized, these events could be a significant source of r- 
process elements. In Ref 9 a different evolutionary path was considered. The point 
of depart is a HS with a deconfined core constituted only by u, d quarks. Then, such 
a core shrinks into a more stable and compact u, d, s configuration in a timescale 
shorter than that of the overlying hadronic material, originating a "quark nova" 
The total energy released in the process may reach values as high as 10^^ erg 
and ~ 10~^ Mq of neutron-rich material may be ejected in the explosion. 

The aforementioned events fail in to explain the GRB phenomenology but could 
shed some light on the provenance of elements heavier than those of the iron-peak. 
High densities and temperatures required to produce these elements are usually 
found in the neutrino-driven wind of type II supernovac , although the fine-tuning 
of the wind parameters necessary to explain the observed abundance pattern is still 
an unsolved issue i^.is 

In the present work an alternative possibility is explored by considering a binary 
system in which one of the components is a " strange dwarf . If this star accretes 
mass what will be its new equilibrium state ? The present investigation indicates 
that the star can either evolve in SD branch by increasing its radius or make a jump 
to the HS (or SS) branch by undergoing a collapse in which the strange core mass 
increases at the expense of the hadronic layer. We argue that there is a critical mass 
(~ 0.24M©) below which a jump to the HS branch is energetically more favorable. 
In this case, the released energy emitted mostly under the form of neutrinos, is 
enough to eject a substantial fraction (or almost completely) of the outer neutron 
rich layers, whose masses are typically of the order of (2 — 5) x Mq. This 
paper is organized as follows: in Section II, strange dwarf models and energetics are 
presented, in Section III, the ejection of the envelope and abundances are discussed 
and, finally, in Section IV the main conclusions are given. 

2. Strange dweirf models 

A sequence of equilibrium (non-rotating and non-magnetic) models were calculated 
by solving numerically the Tolman-Oppenheimer-Volkoff equations^^'^^ (G = c = 



1), e.g. 



dp 
dr 



\p{r) + e(r)][m(r) + 47rr^p(r)] 
r(r — 2m(r)) 



(1) 



and 




(2) 
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The deconfined core is described by the well known MIT bag model from 
which one obtains respectively for the pressure and energy density 

and 

^ = ^ + JLif^Mf^} - Imj) - Imjlni^l^)] , (4) 

where B is the bag constant, here taken to be equal to 60 MeVfm^'^, kf is the Fermi 
momentum of particles of mass m/ and At/ = ^fc j + m j . The sum is performed 
over the flavors f = u, d, s, whose masses were taken respectively equal to m„ = 5 
MeV, TOd = 7 MeV and rus = 150 MeV. 

The hadronic layer begins when the pressure at the core radius reach the value 
corresponding to the density pd of the neutron drip. The equation of state used for 
this region is that calculated in Ref. 25. Notice that the bottom of hadronic layer 
does not represent a true phase transition since Gibbs criteria arc not satisfied. The 
strange matter core may absorb hadrons of the upper layers if they get in contact, 
which is precluded by the strong electric field, as already mentioned. The overall 
equation of state is shown in Fig 1. At this point, it is important to emphasize the 
following point. The transition between the two phases (deconfined and hadronic) 
occurs when the Gibbs conditions (equality between the chemical potential and 
pressure of both phases) are satisfied, case of a first order phase transition. A mixed 
phase has also been proposed but the allowance for the local surface and Coulomb 
energies may render this possibility energetically less favorable In the literature, 
hybrid stars are those with a deconfined core whose transition to the hadronic crust 
is of first order or have a mixed phase. In the present context, we also call hybrid 
stars compact configurations (radii of few km) with a quark core and an outer 
hadronic layer, separated by a strong electric field, as in the case of strange dwarfs, 
but this is clearly an abuse of language. 

The sequence of strange dwarfs and hybrid models was calculated by varying 
the central energy density. In Fig 2 we show the derived mass-radius (M-R) relation 
for our models. The solid curve represents the SD branch, whereas the dashed 
curves represent respectively the branches of pure white dwarfs (on the right) and 
of strange & hybrid stars (on the left). 

The point C in the M-R diagram corresponds to the stability edge characterized 
by a mass M=O.8OM0 and a radius R=1397 km. The point A indicates the position 
of the minimum mass of a stable strange dwarf and coincides with the minimum 
mass model for hybrid configurations. It corresponds to a mass of 0.0237 M0 and 
a radius of 341.82 km. As the central density is increased, the star moves along 
the segment A — > B in the M-R diagram. In this range, strange dwarfs have a 
gravitational mass slightly higher than that of hybrid stars of same baryonic number, 
allowing the possibility for a transition from the SD branch to the HS branch. This 
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Fig. 1. The adopted equation of state describing the deconfined core and the hadronic crust. 
Both regions are connected at the neutron drip density. 

is not the case for strange dwarfs above point B, corresponding to a mass of 0.23 
Mq , since their gravitational masses are smaller than those of HS stars having the 
same baryonic number. The position of the fiducial points A, B and C in the M-R 
diagram as well as the curve AC itself depend on the adopted value for the bag 
constant. A higher value would reduce the mass and radius corresponding to point 
A and similarly, points B and C would be displaced toward smaller masses. This 
occurs because the role of the strong forces increase, leading to more compact core 
configurations ^. 

Physical properties of some SD and HS models are given in Table I. Models in 
both branches are characterized by a given baryonic mass, shown in the first column. 
The gravitational mass (in solar unit) and radius (in km) for strange dwarfs are 
given respectively in the second and third columns, whereas the same parameters 
for hybrid stars arc given in columns four and five. The last column of Table I, gives 
the energy difference AE — (Mj^ — Mq^)c^ between both branches. It is worth 
mentioning that AE is the maximum amount of energy which could be released 
in the process. The variation of the gravitational energy is higher but it covers 
essentially the cost of the hadronic matter conversion onto strange quark matter. 
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Fig. 2. Mass-radius diagram for strange dwarfs (solid curve). Other branches (dashed curves) 
correspond to white dwarfs (WD), strange stars (SS) and hybrid stars (HS). 



Notice that A£ > for masses lower than ~ O.23M0 and A^J < for masses higher 
than the considered limit, as mentioned above. The maximum energy difference 
occurs around ~ O.ISM©, corresponding to Ai? ~ 2.9 x 10^° erg. Strange dwarfs 
above point B in the M-R diagram, if they accrete mass, will evolve along the 
segment B — > C, decreasing slightly the mass and the radius of the deconfined core, 
but increasing slightly the extension of the hadronic layer. The core properties for 
the same models are shown in Table II. Inspection of this table indicates that SD 
along the segment A — > C in the M-R diagram have slightly decreasing deconfined 
core masses and radii. On the contrary, in the HS branch, the deconfined core 
develops more and more as the stellar mass increases. 

A comparison between energy density profiles for SD, HS and SS configurations 
is shown in Fig. 3. All stars have the same baryonic mass [Mb = O.IOGQGA/q). 
It is also interesting to compare our results with those of Rof 6, who have also 
performed similar calculations, but with a slight different equation of state for the 
quark matter. For their model sequence using the same bag constant and the same 
density for the core-envelope transition point, they have obtained comparable values 
for the fiducial points defining the SD branch, e.g., a mass of 0.017 Mq and a radius 
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Fig. 3. Energy density distribution for a strange dwarf (solid line), hybrid star (dotted line) and a 
pure strange star (dashed line) . All configurations have the same baryonic mass, Mb = 0. 1O696M0 . 



Table 1. Properties of Strange Dwarfs and Hybrid Stars. The last 
model corresponds to the minimum mass star and, consequently, has 
only one possible configuration. 



Mb/Mq 




RsD 




Rhs 


AE 






(km) 




(km) 


(xl0-''"crg) 


0.40022 


0.36407 


3582 


0.36577 


8.873 


-30.5 


0.30543 


0.27782 


3547 


0.27878 


8.105 


-17.2 


0.25469 


0.23165 


3505 


0.23179 


7.721 


-2.50 


0.20258 


0.18423 


3442 


0.18411 


7.396 


-1-2.22 


0.16808 


0.15282 


3357 


0.15265 


7.211 


-1-2.93 


0.10696 


0.09723 


3119 


0.09711 


7.249 


-f2.04 


0.05185 


0.04709 


2556 


0.04699 


9.324 


-1-1.74 


0.03626 


0.03291 


2309 


0.03285 


15.02 


-1-1.06 


0.03000 


0.02718 


1756 


0.02718 


28.77 


+0.02 


0.02613 


0.02367 


341.8 


0.02367 


341.8 






of 450 km for point A and a mass of 0.96 Mq and a radius of 2400 km for point 

C. These differences, on the average 20% in the mass and 30% in the radius, are 
probably due to differences in the treatment of the quark matter since the equation 
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Table 2. Core properties of hybrid stars 



Mb/Mq 


M^-/Mq 


R^"-- (km) 




R-S- (km) 


0.40022 


0.01972 


2.663 


0.36560 


8.288 


0.30543 


0.02036 


2.711 


0.25832 


7.135 


0.25469 


0.02117 


2.750 


0.23167 


6.471 


0.20258 


0.02260 


2.801 


0.18401 


5.923 


0.16808 


0.02270 


2.825 


0.15257 


5.202 


0.10696 


0.02277 


2.842 


0.09675 


4.543 


0.05185 


0.02290 


2.844 


0.04676 


3.599 


0.03626 


0.02296 


2.849 


0.03187 


3.169 


0.03000 


0.02298 


2.851 


0.02687 


3.006 


0.02613 


0.02323 


2.868 







of state for the hadronic crust was taken from the same source. More difficult to 

understand arc differences in the radii of some configurations which may attain a 
factor of three. For instance, their model with a mass of 0.0972 Mq has a radius 
of 10800 km while our calculations for a model of similar mass (0.0973 Mq) give a 
radius of only 3130 km. 

The analysis of the energy budget of the SD and HS branches suggests that 
strange dwarfs in the mass range 0.024 - 0.23 M0 and in a state of accretion may 
jump to the HS branch, releasing an important amount of energy 10^° erg). 
The following step consists to study the energetics between the HS and the SS 
branches. A comparison between models in both branches indicates that only for 
masses less than ^ O.I2M0 the transition HS — > SS is possible. Characteristics of 
some computed strange star models are given in Table III. The first three columns 
give respectively the baryonic, the gravitational masses and radii. The last column 
gives the energy difference AE = {Mq^ — Mq^)(?. 



Tabic 3. Strange stars properties. 



Mb/Mq 


m^I'='/Mq 


R*'^' (km) 


AE (xlO'^^erg) 


0.40022 


0.39087 


7.024 


-451 


0.30543 


0.27672 


6.424 


-320 


0.25469 


0.23543 


6.058 


-67.3 


0.20258 


0.18521 


5.608 


-21.4 


0.16808 


0.15278 


5.305 


-0.57 


0.10696 


0.09707 


4.552 


-1-0.80 


0.05185 


0.04699 


3.609 


-1-1.84 


0.03626 


0.03284 


3.207 


-M.23 


0.03000 


0.02711 


3.031 


-1-1.18 


0.02613 


0.02367 


2.876 


-1-0.15 



3. The neutron-rich envelope 
3.1. The ejection mechanism 

The astrophysical environment in which slow and rapid neutron capture reactions 
take place is still a matter of debate. The ejecta of type II supernovae and binary 
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neutron star mergers are possible sites in which favorable conditions may develop. 
Difficulties with the electron fraction in the neutrino-driven ejecta were recently 
reviewed in Ref. 21. For instance, a high neutron-to-sccd ratio, required for a suc- 
cessful r-process, is obtained if the leptonic fraction is small, a condition not 
generally met in the supernova envelope. Non orthodox issues based on neutrino os- 
cillations between active and sterile species, able to decrease Y^, have been explored 
and here another alternative scenario is examined. 

It is supposed a binary system in which one of the components is a strange 
dwarf. The evolutionary path leading to such a configuration does not concern the 
present work. As it was shown in the previous section, strange dwarfs with masses 
in the range 0.024 < M/Mq < 0.24, in a state of accretion, may jump to the HS 
branch since this transition is energetically favorable. These considerations are based 
on binding energies calculated for equilibrium configurations and future dynamical 
models are necessary to investigate in more detail this possibility. 

The jump SD HS is likely to occur within the free-fall timescale, e.g., td ~ 
l/\/Up, which is of the order of a fraction of millisecond. During the transition, 
hadronic matter is converted onto strange quark matter. This conversion leads to 
an important neutrino emission, via weak interaction reactions, consisting the bulk 
of the energy released in the process, which amounts to about 3 x 10^° erg. The 
typical energy of the emitted neutrino pair is about 15-17 MeV, corresponding 
approximately to the difference between the energy per particle of the ud and the 
uds quark plasma. These neutrinos diffuse out the core in a timescale of the order 
of 0.1s (see, for instance, a discussion in Ref 10) through the remaining hadronic 
layers, placed above the high voltage gap, providing a mechanism able to eject the 
outer parts of the envelope. Masses of the hadronic crust are around 2 x 1O~"^M0 
and their ejection requires a minimum energy of about 8 x 10^^ erg, corresponding 
to about 4% of the available energy. 

Neutrinos interact with the crust material through different processes: scattering 
by electrons and nucleons and capture by nucleons. Cross sections for these different 
interactions can be found, for instance, in Ref 22. The dominant process in the crust 
is by far the neutrino-nucleon scattering, whose cross section is 

._=4xlO--Ar^(^)^cm^ (5) 

where E,^ is the neutrino energy and N is the number of neutrons in the nucleus. As 
we shall see below, nuclei in the crust have typically A'' 35 and A ~ 60. Therefore, 
the "optical depth" for neutrinos is 

T^ = J a^-n{j^)ds = 3.3 X 10"^^ y p ds (6) 

The outer hadronic layers have column densities typically of the order of 
(1.3 - 3.0) X 10^^gcm~'^, leading to optical depths of the order of Tj, 0.043- 
0.10, corresponding to a fraction of scattered neutrinos of about 4.2-9.5 %. Thus, 
the momentum imparted to nuclei is able to transfer enough energy to expel the 
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envelope. However, a firm conclusion must be based on a detailed analysis of the 
momentum transfer by neutrinos, coupled to hydrodynamic calculations. 

3.2. Ejected abundances 

The equation of state and the chemical composition of the external hadronic matter 
for densities below the neutron drip were calculated by different authors 24,25,26 _ 
clei present in the hadronic crust are stabilized against /3-decay by the filled electron 
Fermi levels, becoming more and more neutron-rich as the matter density increases. 
The dominant nuclide present at a given density is calculated by minimizing the 
total energy density, including terms due to the lattice energy of nuclei, the energy 
of isolated nuclei and the contribution of degenerate electrons, with respect to the 
atomic number Z and the nucleon number A. 

For a given model, once the crust structure is calculated from the equilibrium 
equations (see Section 2), the mass under the form of a given nuclide (Z,A) can be 
calculated from 



and the integral limits correspond to the density (or pressure) range where the 
considered nuclide is dominant. These nuclides and their respective density range 
were taken from tables given by Refs 25 and 26. Both set of computations have 
used similar mass formulas but slightly different energy minimization procedures. 
As a consequence, some differences in the abundance pattern can be noticed. In 
particular, is the dominant nuclide at densities ~ OApd according to Ref 25, 

whereas in the calculations by Ref 26 the dominant nuclide is 4oZr^'^^. 

When the envelope is ejected, the neutron-rich nuclei are no more stabilized 
and decay into more stable configurations. Notice that the cross section ratio be- 
tween neutrino capture and scattering is ~ Ca/o's ~ 0.008, indicating that neutrinos 
will not affect significantly the original abundance pattern. Nuclei stability were in- 
vestigated using a modified Bethe-Weizsacker mass formula given in Ref 19, more 
adequate for neutron-rich nuclei, and nuclide tables given in Ref 20. 

The resulting masses in the crust for different nuclides are given in Table IV 
and Table V, corresponding to the dominant nuclide data by Ref 25 and Ref 26 
respectively. For both cases, the envelope mass is 3.6 x 1Q~'^Mq. In the first column 
are given the nuclides present in the crust at high pressures, stabilized against 
/3-decay by the presence of the degerate electron sea. In the second column are 
given the stable nuclides originated from the decay of the unstable neutron-rich 
nuclides. The corresponding masses in the envelope are given in the third column 
and abundances by number relative to 2&Fe^^ are given in the fourth column. The 
last column gives an indication of the expected origin of these (stable) nuclides 
in nature: s- and/or r- process and SE for stellar evolution processes in general, 
including explosive nucleosynthesis. 
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Initial 


Final 


Meject 


X/Fe 


origin 




2eFe^^ 


58 


1.000 


SE 




2,Fe^' 


3 


0.050 


SE 


28Ni^^ 


28iVi62 


96 


1.495 


SE 




28iVi64 


55 


0.829 


SE 






17 


0.205 


(s,r) 






25 


0.294 


r 






36 


0.413 


(s,r) 






4 


0.033 


(s,r) 


C„120 
38 '-'^^ 




4 


0.032 


(s,r) 


7„122 

40 /sr 




13 


0.103 


r 



Initial 


Final 


Meject 


X/Fc 


origin 


26-F^e56 


26-Fe56 


61 


1.000 


SE 


28iVi«^ 


28iVi64 


151 


2.166 


SE 




345e80 


17 


0.195 


(s,r) 


32Ge82 


345e82 


25 


0.280 


r 




36i^r«4 


36 


0.393 


(s,r) 


36^rii« 




8 


0.062 


(s,r) 


385ri20 




13 


0.099 


(s,r) 






10 


0.075 


r 



Inspection of table IV reveals a peak aroTind nuclides in the mass range 56-64 
(Fe-Ni peak) also found in ejecta of type la supernovae However, the contribution 
to the iron yield in the Galaxy by one of these events is about 10^ times less than a 
single type la supernova. Nevertheless, in spite of the the small mass of the ejected 
envelope, these events could contribute to the chemical yields of some nuclides like 
Se, Kr and Sn, which are usually supposed to be originated from s and r processes. 
Here their origin is completely diverse, since they are the result of the decay of 
neutron-rich nuclides stabilized by a degenerate electron sea present in the hybrid 
star. 

The required frequency of these events, in order that they could contribute 
significantly to the chemical yield of the Galaxy, can be estimated by using the 
procedure by Ref 28. Assuming that all iron in the Galaxy was produced essentially 
by type la supernovae and adopting for Se, Kr and Sn, nuclides which are here 
supposed to be producd by the collapse of a SD, the present abundances given by 
Ref 29, then the required frequency of these events in the Galaxy is about one each 
1500 yr. 
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4. Conclusions 

Gravitational masses for a sequence of models in the strange dwarf, hybrid and 
strange star branches were computed. Results of these calculations indicate that 
there is a critical mass in the strange dwarf branch, M = 0.24 M©, below which a 
configuration of same baryonic number in the hybrid branch has a smaller energy, 
allowing a transition between both branches. 

If a transition occurs, the envelope radius shrinks typically from a dimension 
of about ~ 3200 km to about ~ 7 km, with conversion of hadronic matter onto 
strange quark matter. In this collapse, the released energy is about 3 x 10^" erg 
carried out essentially by UePg pairs with energies typically of the order of 15-17 
MeV. This value corresponds approximately to the energy per particle difference 
between ud and uds quark matter. Our estimates indicate that neutrino-nucleon 
scattering can transfer about 4-9 % of the released energy to nucleons, which is 
enough to expel partially or completely the hadronic crust, having masses typically 
of about of (2 - 5) X lO^^ Mq. 

The ejecta of these events is rich in nuclides of high mass number and could be 
the major source for the chemical yields of elements like Se, Kr, Sn, if the frequency 
of these events in the Galaxy is about one per 1500 yr. 
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